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Abstract 

The present paper reviews recent attempts to study the development of 
glassy behavior in thin polymer films by means of Monte Carlo simulations. 
The simulations employ a version of the bond-fluctuation lattice model, in 
which the glass transition is driven by the competition between an increase of 
the local volume requirement of a bond, caused by a stiffening of the polymer 
backbone, and the dense packing of the chains in the melt. The melt is 
geometrically confined between two impenetrable walls separated by distances 
that range from once to about fifteen times the bulk radius of gyration. The 
confinement influences static and dynamic properties of the films: Chains 
close to the walls preferentially orient parallel to it. This orientation tendency 
propagates through the film and leads to a layer structure at low temperatures 
and small thicknesses. The layer structure strongly suppresses out-of-plane 
reorientations of the chains. In-plane reorientations have to take place in a 
high density environment which gives rise to an increase of the corresponding 
relaxation times. On the other hand, local density fluctuations are enhanced 
if the film thickness and the temperature decrease. This implies a reduction 
of the glass transition temperature with decreasing film thickness. 
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I. INTRODUCTION AND OVERVIEW 



Many recent studies deal with the influence of confinement on the dynamic behavior of 
glass forming liquids Besides the practical importance of such systems (polymer films 
as protective coatings, flow through porous materials, etc.) one motivation for this research 
is that it might also provide a better understanding of the glass transition in the bulk 
HUH!]- If the glass transition was driven by an underlying correlation length £ which grows 
with progressive supercooling towards the transition temperature T g , confinement should 
truncate the growth as soon as £ becomes comparable to the size, L, of the restrictive 
geometry. Marked deviations from the bulk behavior are then expected. 

However, the trend of these deviations is not obvious. Should the dynamics be acceler- 
ated or slowed down? Theoretical arguments may be put forward for both scenarios. For 
instance, one could speculate that the glass transition is caused by some subtle kind of static 
order which characterizes the spatial arrangement of the particles in the glassy phase and 
spontaneously develops during supercooling. The correlation length should then measure 
the growth of these (solid-like) regions. If this was true, one could invoke an analogy to the 
theory of spin glasses and suggest that the structural relaxation time r scales as r ~ f 
(z being a dynamic critical exponent). When £ ~ L, r should level off and the dynamics 
should become faster in confined geometry. On the other hand, one can also adopt the 
point of view that supercooling makes the dynamics of a glass former become progressively 
cooperative. Cooperativity means that the displacement of a particle is predicated upon 
the simultaneous rearrangement of many neighboring particles and the correlation length 
measures the spatial extent for these cooperative processes. There is evidence for such a 
dynamic heterogeneity from recent simulations anc ^ experiments [ 13 1G ] . In spatial 

confinement some of these processes should be suppressed. Therefore, one expects r to 
increase and the dynamics to slow down compared to the unrestricted bulk. 

A variety of experiments PJT7H2lf1 , computer simulations P,^9|-|37|| and theoretical ap- 
proaches [Q,[3^-|4(J have attempted to display the phenomenology and to elucidate the un- 
derlying mechanisms of dynamics in confinement. The systems studied range from simple 
liquids [ PU| , PT| - P^ , |ID[ ] , over molecular and hydrogen-bonded liquids [|IT|- pT] , p5| ] to polymers 
?8t p6| , p7tp9| . The geometries considered involve three-dimensional cavities pT| , |3T| , |38| , 



pores [pT|- [2l]|P^J55[] , nano-sized fillers embedded in polymer melts [07]], and thin films 
2^- p0tp2| , p3tp6| , p9t^0| , exhibiting different interactions between the glass former and the 



walls of the confinement. These interactions and the arrangement of the molecules close to 
the walls introduce boundary effects which perturb the bulk structure. These effects can 
reinforce or mask the influence of pure confinement on the dynamics. Therefore, a good 
knowledge and control of the surface properties is important for the interpretation of the 
results. 

Such knowledge and control is typically available in computer simulations. Simulations 
study precisely defined - though often highly idealized - systems to reveal possible conse- 
quences of various kinds of confinement. A special kind of confinement which may only be 
realized in computer simulations is a size variation of the simulation box while maintaining 
periodic boundary conditions. This means the following: Usually, the simulated system is 
contained in a cubic box which is replicated in all spatial dimensions. If a particle leaves the 
box on one side, an identical image particle simultaneously enters the box from the opposite 
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side. The simulation cell can thus be thought of as a section of a macroscopic system. The 
fact that the adjacent boxes are identical copies does not influence the physical properties of 
the system if the box size is large (and if there are no such complications, such as long-range 
interactions, spontaneous ordering close to phase transitions, etc.). However, if it shrinks, 
finite-size effects may occur. The development of such finite-size effects during supercool- 
ing has been observed in simulations of binary Lennard- Jones mixtures ||41|| , of hard sphere 
mixtures pZ| , fl3"|] and of silica P^ji5| , and was also treated analytically [B8|. Generally, one 



finds the dynamics to be slowed down with shrinking box size (see |4T| for a comparative 
discussion of the results from Refs. [[Hfl and P2fl). 

An important property of these studies is that the confinement is weak. There is no 
specific interaction with the boundary and no spatial arrangement close to the surface of 
the simulation box. Changes of the bulk structure are very mild or not present at all. The 
same situation is also realized in recent simulations of the relaxation dynamics in spherical 
cavities [jST] and pores p3) , pS|| . In these studies, a sphere or a cylinder is cut out of a large 
bulk sample. Subsequently, only the inner particles are allowed to move, whereas the outer 
particles are frozen. They realize the walls of the confinement. Since the walls are rough 
and adapted to the liquid structure, a fluid particle can be trapped in cavities, into which 
it perfectly fits. Relaxation out of these traps is thus sluggish so that one finds a slowing 
down of the dynamics close to the walls. This slow relaxation also retards the dynamics of 
particles in the inner part of the liquid [j34,46 . 



While the impact of confinement on the structure is negligible in these studies, wall- 
induced perturbations of the bulk structure occur in simulations of thin film geometries. If 
the film geometry is realized by two free surfaces (freely standing films), the density of the 



liquid decreases from the bulk value to zero across the liquid/ vacuum interface [p9| , |30| , ^7 



The interface sharpens with decreasing temperature. On the other hand, if the glass former 
is embedded between two impenetrable and smooth walls, there are pronounced density 
oscillations starting with a value larger than the bulk density at the wall. These oscillations 
can decay towards the bulk density with increasing distance from the wall or propagate 
through the film, depending on the film thickness and temperature considered. These struc- 
tural differences also affect the dynamic behavior. For the freely standing films one generally 
finds that the particles closer to the interface are more mobile than inner particles which 
behaves bulk-like [p9| , |30| , pETf . Therefore, the overall dynamics of the film is faster than the 
bulk at the same temperature. This should lead to a decrease of the average glass transition 
temperature. The same interpretation is also suggested by an analysis of recent experimental 
results for freely standing polystyrene films [^6| and by comparitive computer simulations 
of freely standing and supported polymer film models |HB . 



Whereas the observed speeding up of the particle motion close to a free interface is 
intuitively expected, the situation is not so clear for glass formers confined between two 
completely smooth, impenetrable walls which exert no preferential attraction on the par- 
ticles, but merely represent a geometric confinement. Simulations of glassy polymer films 
5^,50] suggests an acceleration of the local relaxation dynamics, whereas simulations of bi- 



nary hard sphere mixtures rather find a slowing down |32| , |33| . A tentative explanation of 
this difference could be as follows: The aforementioned density oscillations are much more 
pronounced for the binary mixtures than for the polymer films. This is a consequence of the 
model parameters used in the simulations. For the mixtures there seem to be densely filled 
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layers of particles stacked on top of each other at low temperature. Motion has to occur in 
this very high density environment which might therefore be slower than in the bulk and in 
the simulations of the polymer films. 

Thin film geometries are also extensively studied experimentally, especially for polysty- 
rene melts ||22|-|28|. The geometries investigated range from freely standing films [p3}-|26|| 
to supported films (two inequivalent surfaces: SiO"2 and air) p2| , ^4| , p5| , ^8| and polystyrene 
embedded between two SiC>2 substrates, realized either by capping the film with another 
silicon-oxide layer [f24| or by intercalating the polymer in layered-silicate hosts | 27fl . A 



systematic variation of film thickness and chain length has revealed many interesting results 
which have been summarized in recent comprehensive reviews [^T,52 and elicited theoretical 
explanations |39| . 



The present paper also deals with the glass transition of a thin polymer film. It presents 
Monte Carlo simulation results for a simple model of a non-entangled polymer melt confined 
between two solid walls. Its main focus is to investigate static and dynamic properties when 
polymer films of different thicknesses are progressively supercooled. The paper is organized 
as follows: Section describes the background of the model and the employed simulation 
technique. Sections |T| and [IV| illustrate the influence of spatial confinement on the structure 
of the melt and the resulting consequences for the dynamic behavior. The last section [V] 
summarizes and discusses the results. 



II. COARSE-GRAINED LATTICE SIMULATIONS OF POLYMER FILMS 

The present approach uses a lattice model: the bond-fluctuation model This 
model is intermediate between a highly flexible continuum treatment and standard lattice 
models of polymers |56|-|58f|. With the latter it has in common the simple lattice structure 



which is very efficient from a computational point of view ||59|| . However, it differs from them 
in that the set of bond vectors is not limited by the coordination number of the lattice, but 
much larger. In this respect, it resembles more a continuum model. 

A monomer of the bond-fluctuation model does not directly correspond to a chemical 
monomer. It should rather be thought of as representing a group of chemical monomers 
(comprising typically 3 to 5 monomers for simple polymers, such as polyethylene [{3?| ,|60|1 ) . 



A lattice bond should thus be interpreted as the vector joining these groups. This coarse- 
grained vector can fluctuate in length and direction to a much larger extent than its chemical 
counterpart. The bond-fluctuation model accounts for this flexibility of the coarse-grained 
bond vector by associating a monomer not with a single lattice site, but with a unit cell of 
a simple cubic lattice. The monomers can be connected by 108 different bond vectors which 
are chosen such that local self-avoidance of the monomers and uncrossability of the bond 
vectors during the simulation are guaranteed. 

Energy Function and Geometric Frustration. In addition to excluded volume inter- 
action and chain connectivity an energy function 7i(b) is introduced for the bond vectors 
b. It favors bonds of length 6 = 3 and directions along the lattice axes {Ti.(b) = 0) in 
comparison to the rest of available bond vectors (H(b) = e) |6lH63|. Figure |l] illustrates the 
effect of this two-level Hamiltonian. If temperature decreases, each bond attempts to adopt 
the ground state 7~C{b) = 0. A bond in the ground state blocks four lattice sites which can 
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no longer be occupied by any other monomers due to the excluded volume interaction. This 
reduces the amount of accessible volume and generates a competition between the internal 
energy of a bond and the local arrangement of other monomers around it. This competition 
forces some bonds to remain in the excited state. They are geometrically frustrated |]62| , |63 



For instance, the bond of the lower chain in Fig. [I] could only reach the ground state if 
the monomer of the upper chain moved away. However, this motion is only possible if the 
constraints acting on this monomer are released which might in turn require the motion of 
further distant monomers. The relaxation of local geometric frustration is thus predicated 
upon the cooperative rearrangement of many monomers, which can become very sluggish 
due to mutual blocking at low temperatures. Therefore, the development of the geometric 
frustration during the cooling process causes the glassy behavior of the model. It is also the 
driving force which the Gibbs-Di Marzio theory makes responsible for the glass transition 
of polymer melts [ 64 — 166 1 . 



Monte Carlo Methods: Real versus Artificial Dynamics. The competition between 
energetic and packing constraints leads to a slow, glass-like structural relaxation of the melt 
if the usual bond-fluctuation dynamics is used. This dynamics consists of the following 
steps: First, a monomer and a lattice direction are chosen at random. Then, it is checked 
whether the targeted lattice sites are empty (excluded volume interaction) and whether the 
new bonds belong to the allowed set (maintenance of chain connectivity). If these condi- 
tions are satisfied, the energy change AE associated with a monomer displacement in the 
chosen lattice direction is calculated. The attempted move is accepted with probability 
min(l,exp [— 13 AE}) (Metropolis criterion J5M57]]), where (3 = is the reciprocal tem- 

perature (the temperature is measured in units of e/k B ). These local moves are supposed to 
mimic a random force exerted on a monomer by its environment. They lead to Rouse-like 
dynamics which is typical of short polymers in dense melts |68|-[70|] . 



The local dynamics gives rise to a strongly protracted structural relaxation of the melt. 
Although this is exactly the physical phenomenon that we are interested in, such a slow 
relaxation is very disadvantageous to quickly equilibrate the system during the cooling pro- 
cess. To circumvent this problem one can exploit the property of the Monte Carlo technique 
that the elementary move may be adapted at will. Since the final equilibrium state is inde- 
pendent of the way by which it was reached, the realistic local dynamics may be replaced 
by an artificial one which uses non-local moves. A non-local move involves many (or even 
all) monomers along the backbone of a chain. 



The so-called "slithering snake dynamics" is an example for such a collective move |56| , |58 
In this dynamic scheme, one tries to attach a bond vector to one of the ends of a polymer. 
Both the vector and the end monomer are randomly chosen. If the attempt does not violate 
the excluded volume restriction, the move is accepted again with probability exp[— AE/k^T]. 
But now, AE represents the energy difference between the newly added bond and the last 
bond of the other end of the chain, which is removed when the attempt is accepted. Whereas 
the local dynamics propagates the chains by one lattice constant for every accepted move, 
the slithering-snake dynamics shifts the whole chain by one lattice constant. Therefore, one 
expects this algorithm to be faster by a factor of the order of the chain length already at 
high temperatures. The simulations show that this expectation is borne out and furthermore 



that the algorithm is very efficient in equilibrating low temperatures 71.72 
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Simulation Parameters. In the present study a chain always consists of N = 10 
monomers. When taking into account that a lattice monomer roughly corresponds to a 
group of three to five chemical monomers (see above) our simulation deals with fairly short, 



non-entangled |59| , |73| , |7^| oligomeric chains. The number of chains, K, in the rectangular 
simulation box (of size Lx Lx D) is adjusted so that the volume fraction of occupied lattice 
sites is cf) = 8NK/ DL 2 = 0.53. This value is a compromise between two requirements: it 



is high enough so that the model realizes the typical properties of dense melts [59 1, and 
low enough to allow a sufficient acceptance rate of chain moves to equilibrate the melt by 
slithering- snake dynamics [7T| . 

The dimension of the simulation box parallel to the walls is L = 60 (in units of the lattice 
constant; see Fig. [[]). Since the maximum end-to-end distance is R max = 27 (= 9 x 3 = 
(N — 1) x bond length in the ground state), we have L/2 > R max so that there are no finite- 
size effects (this has been tested explicitly by simulations with L = 120). Impenetrable, 
completely smooth and structureless walls are situated at z = 1 and z = D. Technically, 
these walls are generated by disallowing monomer moves to z < 1 or z > D. The film 
thickness ranges from D = 6 (~ 1.5_R g ; R g = bulk radius of gyration) to D = 60 (~ 15i? g ). 

In order to improve the statistics several independent simulation boxes are treated in 
parallel. The total statistical effort involves between 43200 and 86400 monomers, depending 
on the system under consideration. This allows us to keep the numerical uncertainties of 
the results at a low level. Unless otherwise stated, the statistical errors are always of the 
order of the size of the symbols in the following figures. 

III. STATIC PROPERTIES: GYRATION TENSOR AND DENSITY PROFILES 

Structural properties represent an important input for analyzing the dynamics of a poly- 
mer melt (and of any other system in general). Close to a solid interface the structure of 



the melt markedly deviates from the behavior of the unconstrained bulk [ff5|,|76fl . This is 



pointed out by analytical approaches |77|-|88l, computer simulations of lattice |89| |95| and 



continuum models [^ , [?^ , [7^ . ^ - |83| , |96Hl02|| , and recent experiments ||103|| (however, see also 



10 1|| for different experimental results). The present section describes the influence of hard 



walls on the static properties of our model. 

A. Gyration tensor 

A polymer chain in the bulk can adopt a multitude of different configurations. Each 
configuration has a certain distribution of monomers around its center of mass. This distri- 
bution determines the instantaneous shape of the chain. Classical mechanics suggests that 
a quantitative measure of this distribution may be provided by the moment of inertia tensor 



p5] 



= (TrQ)l-Q , 

where 1 denotes a 3 x 3 unit matrix and Tr Q is the trace of the tensorial generalization of 
the radius of gyration (note that (Tr Q) = R 2 ) 
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1 N 

Qaf3 = J7 ^2 [ ( r "> Q ~ R cm,a) (?"n,/3 ~ Rcm,p) («, /? = 1, 2, 3) . (1) 

n=l 

Here, r n Q , and R cra ,a are the ath spatial component of the position vectors to the monomer 
n (r n ) and to the center of mass (i? cm ), respectively. 

The average of or equivalently of the gyration tensor Q are good indicators of the 
characteristic shape of a polymer because they reflect the average distribution of monomers 
in the internal coordinate system of a chain. This coordinate system is given by the principal 
axes obtained after diagonalization of Q. The corresponding eigenvalues, Ai, A2 and A3, 
measure the length of the axes. If the (instantaneous) monomer distribution was spherical, 
all eigenvalues would be equal. 

However, the seminal work of Sole and Stockmayer already pointed out that the average 
shape of a (single isolated) polymer is far from spherical ||106|| . This conclusion has been 
corroborated by various further studies ||105J107| - tL09|| . All principal axes of a polymer have 



different lengths. For a random walk the ratio of the average eigenvalues is given by: (Ai) : 



(A 2 ) : (A 3 ) = 12.07 : 2.72 : 1 ||109|| . Therefore, even a random walk is distorted with respect 
to a sphere: Its size is extended along the largest principal axis and reduced in directions 
of the axes corresponding to (A2) and (A3). Since a random walk is a viable model for (the 
large scale properties of) a chain in the melt, the average shape of chains in the melt rather 
resembles a flattened ellipsoid than a sphereQ. 

In addition to the eigenvalues two other quantities were introduced to discuss deviations 
from a spherical structure: the asphericity A and the prolateness S ||105| , |110|| . They are 
defined by 

3 _ 2 

3 TrQ 2 1 

and 

detQ n(Aa-A) 

5o = 27^*3=^^ (3) 

(TrQ) 3 A 1 ' 

with A = (Ai + A 2 + A 3 )/3. 

For a sphere one has A a = A so that A = So = 0. Two extreme deviations from this high 
symmetry can be considered: On the one hand, the sphere may be squashed to a disk. Then, 
Ai = and A 2 = A 3 = 3 A/2 which yields A = 1/4 and S = —1/4. On the other hand, it 



1 Recent studies | 10§| , 109 of the shape of Gaussian chains indicate that the visualisation of a 
random-walk polymer as a flattened ellipsoid is not completely correct. The density distribution of 
monomers in the coordinate system of the principal axes exhibits a slight minimum at the origin 
(i.e., at the center of mass) for the largest axis, whereas it has a maximum for the other two axes. 
Therefore, the shape is rather dumbbell-like. 
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can be stretched to a rod. Then, Ai = 3A, whereas A2 = A3 = so that Ao = 1 and So = 2. 
Therefore, the asphericity only measures deviations from the spherical structure, whereas 
the prolateness additionally determines by its sign whether the deviation is disk-like, i.e., 
"oblate" (S < 0), or elongated, i.e., "prolate" (S > 0). 



B. Asphericity and Prolateness of Polymer Films 

High- Temperature Results. Figures ^| and [3] compare the variation of (Ao(z cm )) and 
(So(z cm )) with the distance of the chain's center of from the (left) wall for different 

film thicknesses D at T = 00. Here, (•) denotes the average over all chains in the system. 
For T = 00 only excluded volume interactions between the monomers and between the 
monomers and the walls are effective. This temperature is therefore representative of the 
high-temperature liquid state of the melt. In this state, the average shape of chains in the 
inner part of the film resembles a prolate ellipsoid. The numerical values for (A (z cm )) 
and (So(z cm )) are close to those of a random walk [|107|| . The bulk-like inner part extends 



from the middle of film to about z cm ~ 2R g = 7.36, where the interfacial region starts. 
Since the melt is confined between two walls, a bulk-like inner part can only be observed 
if D > 18 (~ 5Rg). For smaller thicknesses (D < 12 « 3R g ) the film just consists of 
interfacial region because the perturbations of the structure, which propagate from both 
walls, interfere. Similar results have already been observed before [ 90 , 91 , |95 , 1 00| , |1 02j , |l 1 1| . 



Let us first consider the larger films D > 18. For these thicknesses the profiles of 
(A (z cm )) and (S (z cm )) are independent of D. When z cm < 2i? g , the chains first slightly 
expand along the longest principal axis, then shrink, pass through a minimum and finally 
become very elongated when their center of mass lies at the wall (i.e., at z = 1). This 
variation of the eigenvalues is accompanied by a reorientation of the principal axes (see [|112|1 ). 
Whereas the chains can orient freely in the bulk-like inner region of the film and thus appear 
spherical on average, the solid wall singles out those configurations, in which the two largest 
axes are aligned parallel to it. The corresponding eigenvalues are bigger than those measured 
in the film center [(Ai(z cm = 1)) ~ 12.81 > (\i(z cm = 15)) ~ 10.23, (A 2 (l)) ~ 2.40 > 
(A2(15)) ~ 2.32], whereas the third eigenvalue is considerably smaller [(As(l)) ~ 0.34 < 
(As(15)) ~ 0.80]. Therefore, the chains are not only oriented parallel to the wall, but also 
distorted. They adopt a rather flattened shape. This influence of the wall on the structure 
of the chains is also confirmed by other simulations piP,pip9H9^,|95|j97|, [T00| , [I02| , [TTT1 . 



A distortion of the structure and orientation of the chains is entropically unfavorable 
because they strongly restrict the number of accessible configurations. As z cm increases, 
the disk-like ellipsoid therefore turns away from the parallel alignment, shrinks in directions 
of the longest axes and expands along the smallest axis. This effect is most pronounced 
for z cm = 3, where the polymer concentration has a maximum (see Fig. |5|). In order to 
accomodate many chains at the same distance each chain has to be compressed so that the 
average shape is least prolate [(Ai(z cm = 3)) ~ 8.86, (A 2 (;z cm = 3)) ~ 2.31, (A 3 (z cm = 3)) ~ 
0.84]. This effect is stronger for D = 6 because the chain density is much higher at z cm = 3 
(middle of this film) than for the thicker film (see Fig. [5]). 

If z cm increases further, the chain expands again and (A (-2 C m)) and (So(z cm )) go through 
weak maximum at z cm = 6 before crossing over to the respective bulk values. The maximum 
corresponds to chains which are on average slightly oriented in direction perpendicular to the 
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wall. Such chains can (presumably) already touch the wall with some of their monomers and 
thus contribute to the large monomer concentration at z = 1 (see Fig. |^). This interpretation 
is corroborated by simulations of liquid n-tridecane, in which the variation of the monomer 
distribution around the center of mass shows exactly this behavior, as z cm approaches the 
wall [ |113|| . Since z cm = 6 = D/2 for D = 12, chains at z cm = 6 can reach both walls, which 
enhances (A (z C m)) and (So(z cm )) for this particular thickness in comparison to the other 
films. 

A similar behavior can also be observed for quantities which are experimentally better 
accessible than (A (,2 cm )) or (S (z cra )). For instance, for the radius of gyration P3||M|] (or the 



end-to-end distance ||114|| ). When measuring the components of radius of gyration parallel 



and perpendicular to the wall, one finds that the parallel component R g ,\\(z cm ) is large at 
-2 cm = 1? whereas R g ,±(z cm = 1) is small (see Fig. f|). With increasing separation from 
the wall the perpendicular component R gt ±{z cm ) increases towards a pronounced maximum 
at z cm = 6 before crossing over to the bulk value. On the other hand, R g: \\(z cm ) has a 
shallow minimum at z cm = 6 only, but otherwise decreases continuously towards the bulk 
value. For both components the bulk value is reach if z cm « 2i? g . This behavior is found 
in many other simulations |79l^pjg9l^p5l,p7i [100| , [l02lplll , can be reproduced by self- 



consistent field theory (for confined binary polymer blends) [ |115[ | and is also suggested by 
some experiments ||103|| (however, [103] reports that R g remains essentially bulk-like, even 
for D/R g « 0.5). 

Dependence on Temperature. The insets of Figs. |] and |3| illustrate the temperature 
dependence of (A (z cm )) and (S (z cm )) for the smallest film thickness D = 6 (~ 1.5R g ). 
Although this film has no bulk-like inner part, the infinite temperature results for the as- 
phericity and pro lateness at z cm = 1 coincide with those obtained for larger D. Therefore, 
the structure of the chains at the wall is independent of film thickness. Other simulations 
support this finding [|90| , |102|1 . When z cm increases, (A (z cm )) and (So(z cm )) decrease and 



become much smaller than the bulk value in the middle of the film. This is (probably) a 
consequence of the high polymer concentration at z cm = 3, as pointed out above (see also 
Fig. §). 

If the film is progressively supercooled, this behavior changes in two respects: First, the 
values of (A Q (z cm )) and (So(z cm )) at the wall increase. This indicates that the chains become 
more prolate with decreasing temperature. Second, the minimum at z cm = 3 turns into a 
maximum of about the same height as found at the walls. Furthermore, there are minima 
between the maxima at z cm = 1,3,5 which have almost the same value as that at the walls 
for T = oo. The resulting zig-zag structure is in phase with the profile of the chain density 
(Fig. [|). Contrary to the behavior at T = oo, chains located at positions of high density are 
extremely prolate. This means that all chains of the film, and not only those at z cm = 1, 
are preferentially oriented parallel to the walls and flattened. 

This temperature dependence of the profiles is a consequence of the model's energy 
function. Remember that this function favors large bond vectors which point along the 
lattice directions and are thus parallel to the walls. At low temperature, the energy function 
will therefore reinforce the influence of the walls to align chains parallel. The combined effect 
should be to single out those chain configurations which almost lie completely in the lattice 
layer next to the wall. Since this layer is rather densely filled (see discussion of Fig. |5|) and 
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thus essentially impenetrable, it represents another (rugged) wall in front of the original one. 
Effectively, the film has become thinner by two lattice constants (one from each wall). This 
reduction of the available volume limits the orientational freedom of the remaining chains 
in the middle of the film so that they also align parallelly. 

C. Density Profiles of Chains and Monomers 

Chain Profiles. The discusssion of the previous section pointed out that the shape of 
chains at the wall is flattened. Such a deviation from the bulk structure and their parallel 
orientation are entropically unfavorable. At T = oo, where only entropic effects are present, 
one can therefore expect the polymers to avoid the vicinity of the walls. Figure [| shows that 
the chain concentration, p c (z cm ), is in fact negligible at z cm = 1 for all film thicknesses. On 
the other hand, if z cm increases, the concentration quickly rises, attains a maximum close 
to z cm = R g and then decreases towards the bulk density which is reach for z cm > 2R g . 
This behavior is found for D > 12 (^ 3i? g ) and confirmed by simulations with other models 
P|,|gTlJg^,|gg|^,PB|,^, |iI?q , |ni^rm[] . The influence of film thickness on p c is weaker than for 
the asphericity and prolateness discussed in the previous section. Large deviations are only 
seen for the thinnest film D = 6, where D is so small that the depletion effects occurring at 
both walls interfere and lead to a high concentration in the middle at z cm = 3. The same 
argument was also put forward to interprete the simulation results of [ |102| . 



If temperature is reduced, the shape of p c {z cra ) alters completely. This is exemplified 
again for D = 6 in the inset of Fig. Since a decrease of temperature favors long bond 
vectors along the lattice axes, the triangle-shaped profile of T = oo with a low concentration 
at the wall gradually turns into an enrichement of chains at z cm = 1. These chains lie almost 
flat in the lattice layer next to the wall. The maximum number of chains which a layer can 
accomodate at low temperatures can be estimated as follows: Let us assume that a chain is 
a rod with iV — 1 = 9 bonds in the ground state (bond length = 3). The end-to-end distance 
is then 27. Due to excluded volume interactions, a 60 x 60 lattice layer can take up 60 chains 
at maximum. So, p™ ax (^ cm = 1) = 0.016. A comparison of this estimate with Fig. [| shows 
that the layer at the wall is very dense. Due to excluded volume interaction the following 
layer is depleted, which in turn allows for a larger concentration in the subsequent layer 
and so on. This zig-zag structure is a result of the interplay between the model's energy 
function, the wall and the underlying lattice. 

Monomer Profiles. The monomer profile p m (z) is defined as the average density of 
monomers which are situated at a distance z from the (left) wallQ. It exhibits an oscil- 
latory structure not only at low, but also at high temperatures (see Fig. |6]). This structure 
can be explained as a result of the competition between packing constraints and loss in en- 
tropy. The loss in entropy is caused by the reduction of accessible chain configurations near 
an impenetrable wall. It produces an effective repulsive force pointing away from the wall. 
This force competes with another effective force exerted by the densely packed chains in the 



2 The position of a monomer is associated with the z- value of the lower left corner of its unit cell. 
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inner part of the film. They tend to push the polymers, which are close to the wall, towards 
the wall. At melt-like densities packing constraints dominate and lead to an enrichment of 
monomers at the wall. Due to the mutual exclusion the monomer concentration is reduced 
in the next layer, which in turn allows an enrichment in the subsequent layer and so on. 
Therefore, the monomer profile is expected to exhibit a sequence of maxima and minima, 
the amplitudes of which should decay to zero when approaching the inner bulk-like portion 
of the film ]^gI|^,gg^^gg- ^ , gTjfIT§ . 



Figure |j] illustrates this behavior of p m (z) at T = oo and T = 0.2 for D = 30 (see also 
At T = oo, the comparison of the low polymer, but high monomer concentration 
at the wall suggests that most of the monomers at z = 1 belong to different chains. This 
conclusion is true, as the calculation of the average number of monomers, which belong to 
the same chain and are in the same layer, shows. The number is about 3 at T = oo, but 
increases to 7 or 8 at T = 0.2 (see also p3| , |112|| ). This is again evidence for the tendency of 



the chains to orient parallel to the wall under the influence of the model's energy function. 

IV. DYNAMIC PROPERTIES OF THE POLYMER FILMS 

Previous work on the dynamics of polymer films focused on the behavior of mean-square 
displacements and related quantities, such as the monomer mobility and the chain's diffu- 
sion coefficient parallel to the walls |3^j9^JTI^| . Another important means to study dynamic 



properties are time-displaced correlation functions. The present section discusses two differ- 
ent kinds of these functions: the incoherent intermediate scattering function which probes 
density fluctuations, and the correlation functions of the Rouse modes which are sensitive 
to reorientations along the backbone of a chain. 

A. Incoherent Intermediate Scattering Function 

Let r m (t) be the position vector to the mth monomer at time t. The incoherent scattering 
function can then be written as 

M 

^ = m E ( ex P fa • M*) - r -(°)]) ) > ( 4 ) 

m=l 

where q = (q x ,q y , 0) denotes a wave vector corresponding to the x, y-plane of the simple 
cubic lattice and (•) stands for both the thermodynamic average and the average over all 
wave vectors with the same modulus q — \q\. Physically, the incoherent scattering function 
measures displacements of a monomer in time. The dominant contribution to the decay of 
<&q{t) comes from motions of the order 27r/g. 

Qualitative Aspects of the Decay. An analysis at high temperature (i.e., at T = oo) 
pointed out that the influence of finite film thickness on the decay of <& s q (t) diminishes with 
decreasing q- value, i.e., when density fluctuations on larger and larger length scales are 
probed p9[. Thus, the present study focuses on the peak of the static structure factor, 



q = 2.94 1 116], which corresponds to the scale of a bond in real space. 
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Figure ^ shows the time-dependence of $ s q (t) at this q- value for various temperatures 
and two different film thicknesses, D = 6 (~ 1.5R g ) and D = 30 (~ 7.5R g ). The figure 
clearly illustrates the slowing down of structural relaxation as the temperature decreases 
from the normal liquid (T = 0.35) to the supercooled state (T = 0.2) of the melt. The term 
"supercooled state" refers to temperatures close to, but still above the critical temperature 



of mode-coupling theory (MCT) P, |117[HjT9f for the bulk, which was estimated as T c « 0.15 

in 

Two observations can be made from Fig. [7]: First, the thin film always relaxes faster 
than the thick film. This difference increases with decreasing temperature. Second, there is 
hardly any influence on the shape of the decay for both film thicknesses as long as T > 0.2. 
If T = 0.2, a shoulder begins to emerge at intermediate times for D = 30. A similar feature 
is not visible for D = 6. This shoulder can be interpreted as the onset of the MCT-/3 process 
whose characteristic signature is a two-step relaxation of & s q (t) [3, |1 17HTT9|| . Physically, 



the /^-process corresponds to the relaxation of particles in "cages" formed by their nearest 
neighbors. As temperature approaches T c , a particle is trapped for some time in its local 
enviroment (= "cage") before it can escape and diffuse to an adjacent cage P, |117HTT9| . The 



approximate theoretical treatment of this picture compares fairly well with simulation data 
for the bulk |6l|j62]rr20|J121|| . The present results suggest that this relaxation behavior might 
also develop in a thin film geometry at lower temperatures than studied. 

Scaling Behavior of the Incoherent Scattering Function. In order to illustrate these 
qualitative properties from a different point of view we determined the relaxation time 
t* from the correlators by posing Q q (r q ) = 0.629 and plotted ® B q (t) versus t/r q . Such a 
representation shifts the curves on top of each other and thus allows a better comparison of 
the influence of temperature and thickness on the shape of the decay. As expected from the 
qualitative discussion above, the data for both film thicknesses collapse at high temperature 
(exemplified by T = 0.35 in Fig. ||). This indicates that the shape is unaffected by the 
confinement. The main influence is a change of the relaxation time. 

However, such a collapse is only possible for the final decay, i.e., the late time a-process, 
if the temperature belongs to the supercooled regime (T < 0.23). The shoulder is now 
clearly visible for T = 0.2 and D = 30 at intermediate times. Exactly in this intermediate 
time window the simulation results for D = 6 markedly deviate from D = 30, whereas 
they remain (very) close to D = 30 if T > 0.2. The data for D = 6 and T = 0.2 do 
not exhibit the shoulder. They rather resemble simulation results of D = 30 obtained at a 
higher temperature than T = 0.2. In other words, the thin film seems to behave like the 
thick film at a larger temperature. Therefore, the present analysis suggests that a possible 
influence of geometric confinement on the dynamics of the polymer films could be a shift of 
the temperature scale to higher values. Evidence in favor of this interpretation is provided 
by simulation results for a bead-spring model of a glassy polymer film |50| and by the 
temperature dependence of the relaxation time T q . 



3 Contrary to the present simulation data, the bulk configurations, from which the estimate T c ~ 
0.15 was derived, were not fully equilibrated. It is possible that a complete removal of these residual 
non-equilibrium effects shifts the result for T c (to slightly higher temperatures). 
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Temperature and Thickness Dependence of the Relaxation Time. Figure ^| shows 
the temperature dependence of the relaxation time r| for the both film thicknesses, D = 6 
and D = 30, and compares it with that of bulk. One can see that the curves gradually 
splay out with decreasing temperature. The films always relax faster than the bulk, and the 
relaxation time of the film is the smaller, the smaller its thickness. If one extrapolates this 
trend to low temperature by a Vogel-Fulcher equation 



s 



TS ,= A >^{f~r ) ' (5) 

one obtains a Vogel-Fulcher temperature Tq which decreases with decreasing film thickness 



49fl . Therefore, density fluctuations on the length scale of bond suggest a reduction of the 



glass transition temperature with decreasing film thickness. 



B. Rouse Modes 

Possible quantities to study reorientations are correlation functions of the Rouse modes 
. The Rouse modes are the cosine transforms of the position vectors, r n , to the mono- 



mers. For the discrete polymer model under consideration they can be written as ||122 

N 

..V ■ y v J 



Xjt) = -y r n (t) cos I " U l( 2)p7r ) , P = 0,...,N-1. (6) 



^X)r n (t) cos(- 

n=l ^ 



One can think of the Rouse modes as standing waves along the backbone of the chain. The 
first mode has nodes only at the ends and is thus sensitive to reorientations of the whole 
polymer. The second mode has an additional node in the middle. It roughly divides the chain 
into two halves and probes reorientations of these two segments. The higher Rouse modes 
further decompose the chain. Approximately, one can consider the pth mode as a quantity 
which is sensitive to orientational dynamics of a chain segment with N/p monomers. For 
the studied chain length N = 10 the fifth Rouse mode thus measures relaxation processes 
on the same length scale of & s q (t) at the maximum of the static structure factor. Therefore, 
we concentrate on this mode. 

The auto-correlation function of the fifth mode is given by 

where r 5 is the relaxation time of the mode. The second equality of Eq. (|7]) is the prediction 
of the Rouse model which is only approximately satisfied for the present model [f74] , |123|l in 
the bulk (and other models as well ||124|| ), if p — 1. For higher mode indices deviations from 



a simple exponential behavior are observed |7^J123|| . So, we expect to find similar differences 



between the simulation data for $5s(t) and exp(— t/r 5 ) for the polymer films studied here. 

Since the Rouse modes are vectors, one can distinguish between a component parallel 
and perpendicular to the wall. In the following the auto-correlation function parallel to 
the wall shall be considered only. Why? Imagine a polymer film with a thickness of one 
monomer. In this two-dimensional geometry no relaxation perpendicular, but only parallel to 
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the walls can occur. In a film of finite thickness one may expect that the parallel component 
still dominates the dynamics of the correlation function (0) provided the thickness does 
not reach bulk-like dimensions. Therefore, we focus our attention on the relaxation of the 
parallel component of the fifth mode, $55 \\(t), in the remainder of this section. 



Qualitative Properties of the Decay. Figure |T0| shows the time-dependence of $5511 (£) 
for the same temperatures and film thicknesses as studied before for the incoherent scattering 
function. Again, several observations can be made: The influence of film thickness on the 
shape of the decay is almost negligible. It is still weaker than for $g(£). Throughout the 
temperature range studied the thick film (D = 30) remains very close to the bulk, and 
$55,11 {t) always decays in a single step for both film thicknesses and for the bulkQ. There is 
no indication of a two-step relaxation contrary to $„(£)■ The film thickness seems to shift 
the relaxation time only. 

This is also illustrated in Fig. |TT] which shows an attempt to scale all data onto a master 
curve by plotting $55,11 (i) versus t/r 5 in analogy to Fig. ||. As in the case of $g(£), the 
relaxation time was defined via $55,11 (T5) = 0.629. This scaling almost perfectly collapses all 
data onto a common curve, the late time decay (i.e., $55,11 (t) < 0.5) being eventually slightly 
different between the supercooled state (T < 0.23) and the normal liquid state (T > 0.35). 
This indicates that the confinement effects visible in Fig. |10| leave the shape of $55,11 (t), but 
only affect the relaxation time. 



Temperature and Thickness Dependence of the Relaxation Time. Figure [L2| shows 
the temperature and thickness dependences of r 5 . Contrary to t*, one finds that r 5 increases 
with decreasing film thickness. The thin film (D = 6) always relaxes more slowly than the 
thick film (D = 30) which is in turn (a bit) slower than the bulk. A tentative explanation 
could perhaps be as follows: The discussion of the static properties showed that chains in 
the film of thickness D = 6 are oriented parallel to the walls. At high temperature they 
crowd in the middle, whereas an oscillatory density profile with a high concentration at the 
walls develops during supercooling. In both cases the decorrelation of X 5 ,|| has to take place 
by many parallel displacements of all monomers of a chain in an environment where the 
chain density is higher than in the thick film or in the bulk. So, the relaxation could be 
slowed down for D = 6 in comparison with D = 30 which is more homogeneous and thus 
more bulk-like. 



V. SUMMARY 

This paper reports simulation results for a simple model of glassy polymer films. The 
model consists of short (non-entangled) monodisperse chains. The monomers of the chains 
solely interact by excluded volume forces with each other and with two completely smooth 
walls that define the thin film geometry. The glassy behavior is brought about by a com- 
petition between the internal energy of a chain, which tries to make the chain expand at 



^This could be a property of the studied model because a two-step decay of the Rouse modes was 



found in a molecular-dynamics simulation of a bead-spring model [ 124 1 
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low temperatures, and the dense arrangement of all monomers in the available volume. The 
temperatures investigated are taken from the temperature interval above the bulk critical 
temperature of mode-coupling theory. 

With this model we studied both static and dynamic features of polymer films of various 
thicknesses. For the interpretation of the static results it is important to realize that the 
instantaneous shape of a polymer resembles an ellipsoid. Close to a hard wall, the ellipsoid 
tends to orient parallel to it (see Fig. |TB| ). If the film thickness is large enough (D > 5R g ), 
there is a bulk-like inner region with free orientation of the chains. But, as the film thickness 
decreases, this orientational freedom becomes more and more limited. If D < 4iL, the 
perturbations of the structure induced by both walls interfere so that chains also align 
parallel to the walls in the inner portion of the film. With decreasing temperature this wall 
effect becomes strongly enhanced as a consequence of the internal energy of the chains which 
favors bond vectors that point along the lattice axes and are thus parallel to the wall. 

The deviations from the isotropic structure of the bulk also affect the dynamic properties 
of the films. The relaxation of vectors, such as the Rouse modes, is dominated by reorien- 
tations parallel to the wall due to the alignment of the chains if the film thickness is small 
or if the temperature is low. Under these conditions (small D and T) reorientations occur 
in an environment of high chain density which might be the reason why the relaxation time 
of orientational correlation functions increases with decreasing film thickness. 

On the other hand, these orientations could also influence local density fluctuations. 
Oriented chains should have a smaller number of contacts with other chains compared to 
the bulk because they are less intermingled. This might be an explanation for the faster 
relaxation of <f) s q {t) in the thin film (D = 6) already at high temperature and particularly at 
low temperature, where orientation effects become more pronounced due to the impact of 
the model's energy function (see also [ |125| | for a similar argument to interprete experimental 
results for freely-standing polymer films). 

At low temperatures, however, another effect could also contribute. Theoretical de- 
velopments [|117| , |119|| as well as tests by experiments and simulations [^. |118| . |T26| - |12^| have 
suggested that a significant contribution to the slowing down of structural relaxation stems 
from the mutual blocking of particles which are close in space (particles in the "cage"). This 
cage effect becomes important for the dynamics in the supercooled state and is determined 
by the liquid structure of the bulk, particularly by the first neighbor shell. Deviations from 
bulk behavior should be observed if this structure is perturbed. The structure is perturbed 
in the polymer films studied. The completely smooth walls of our model cut off the liquid 
structure and thereby remove part of the obstacles (i.e., other particles) which impede the 
displacement of a tagged particle (see Fig. [TJ]). The walls act like a lubricant with respect 
to the bulk, which enhances the mobility of the nearby monomers with respect to the bulk. 
These more mobile monomers can transfer part of this impetus towards the inner portion 
of the film so that the effect should propagate over a certain distance away from the wall 
before it is gradually damped out. 

This effect could complement the orient at ion- induced reduction of contacts for thin films 
in the supercooled state where the dynamics is very cooperative. Even if the chains densely 
populate the layer next to the wall at low temperatures for small film thickness and thus 
essentially form a more or less rugged surface in front of the wall, this surface is not frozen 
in, as assumed in the right panel of Fig. [TJ], for instance. On the one hand, this surface 
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makes adjacent chains to orient parallel to it, thereby effectively reducing the number of 
contacts, and on the other hand, it should be lubricated by the smooth underlying wall and 
pass part of this enhanced mobility over to neighboring layers. Both effects could therefore 
reinforce one another. 

However, there is still another effect which might eventually outweigh this enhanced 
mobility. The discussion of the static properties showed that the monomer density is larger 
at the wall than in the bulk and decays towards the bulk density in an oscillatory fashion 
with increasing distance from the wall. If the external control parameters (temperature, 
pressure, film thickness) are such that the density oscillations become very pronounced, the 
film could be divided into highly populated layers separated by sharp depletion zones. This 
seems to be the case in the above cited simulations of hard sphere mixtures [f32i [33f . At low 
temperature the particle density at the wall is about a factor of 4 larger than the bulk value 
in these studies. In the present simulation, it is only a factor of about 2. It is plausible that 
the mobility should decrease if the density surmounts a certain threshold, which could in 
turn slow down the dynamics of the whole film. Therefore, even for a completely smooth 
wall the extent of the structuration of the liquid close to it might lead to faster or slower 
dynamics compared to the bulk. Qualitatively, this argument is similar to those presented in 
a discussion of the influence of wetting properties on diffusion in confined liquids 



Finally, the preceding discussion always assumed that there is no preferential attraction 
between the monomers and the walls. If such an attraction, however, exists and if it is 
strong, one would expect a rather immobile interfacial layer of adsorbed particles to form. 
For small film thicknesses this layer should considerably influence of the overall film dynamics 
and could contribute to slowing down the structural relaxation of the film. In this case, the 
glass transition temperature should increase for thin films with respect to the bulk. There is 
evidence for this behavior from simulations of both freely standing and supported polymer 



films and from experiments of polymer films [Pq , |131|| (see also [p0| , |132|| for comparable 
results of molecular glass formers). 
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FIGURES 



FIG. 1. Sketch of the simulation geometry (left panel) and of the model (right panel). The 
simulation box is confined by two hard walls in the z-direction, which are a distance D apart 
(D = 6, . . . , 60 « (1.5, . . . , 15)-R g ; R g : bulk radius of gyration). In the x- and y-directions periodic 
boundary conditions are used (exemplified by the bond leaving the bottom and reentering at the 
top). The linear dimension in these directions is L = 60. The right panel shows a possible config- 
uration of two different chains. All bonds have energy e (e/fce = 1: this defines the temperature 
scale) except the bond (3,0,0) which is in the ground state (two-level system). This vector blocks 
four lattice sites (marked by o) due to the excluded volume interaction. This interaction also 
forbids the jump in direction of the arrow. 

FIG. 2. Asphericity (Ao(z cm )) of a chain versus the distance, z cm , of the chain's center of mass 
from the (left) wall at T = oo (main figure). (Ao(z cm )) is defined by Eq. @. Since the profile of 
(Ao(-Zcm)) is symmetric around the middle of the film, only one half is shown. Various film thick- 
nesses D, ranging from D = 6 « 1.5R g to D = 60 w 15i? g , are compared (R g {T = oo) ~ 3.68 = 
bulk radius of gyration; vertical dashed line). The result for the bulk, (Ao) = 0.411, is also depicted 
(horizontal filled circles •). This value is close to that expected for a random walk (Ao) = 0.396 
| p.07[ |. The inset illustrates the temperature dependence of (Ao(z C m)) for the thinnest film D = 6. 
Here, the profile over the whole film is shown. The temperature varies from the high-T, liquid 
state (T = oo) to the supercooled state of the melt (T = 0.18). Lines are guides to the eye only. 
All lengths are measured in units of the lattice constant, and temperature is measured in units of 
e/k B (see Sect. and Fig. |l|). 

FIG. 3. Prolateness (So(z cm )} of a chain versus the distance, z cm , of the chain's center of mass 
from the (left) wall at T = oo (main figure). (So(z cm )) is defined by Eq. (|3|). Since the profile of 
(So(z cm )) is symmetric around the middle of the film, only one half is shown. Various film thick- 
nesses D, ranging from D = 6 as 1.5i? g to D = 60 « 15i? g , are compared (R S (T = oo) ~ 3.68 = 
bulk radius of gyration; vertical dashed line). The result for the bulk, (So) = 0.499, is also depicted 
(horizontal filled circles •). This value is close to that expected for a random walk (So) = 0.481 
|]107|| . The inset illustrates the temperature dependence of (So(z cm )) for the thinnest film D = 6. 
Here, the profile over the whole film is shown. The temperature varies from the high-T, liquid 
state (T = oo) to the supercooled state of the melt (T = 0.18). Lines are guides to the eye only. 
All lengths are measured in units of the lattice constant, and temperature is measured in units of 
e/k B (see Sect. and Fig. 0). 
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FIG. 4. Profile of the radius of gyration, R 2 (z cm ), measured parallel and perpendicular to the 
wall, at T = oo (main figure). z cm denotes the distance of the chains' center of mass from the left 
wall. Since the profiles are symmetric around the middle of the film, the left half is only shown. 
Three different thicknesses are presented: D = 6 (~ 1.5R g ), D = 12, and D = 30 (~ 7.5R g ). 
The bulk radius of gyration, R g (~ 3.68), is indicated as a vertical dashed line. All lengths are 
measured in units of the lattice constant. The inset shows the counterpart of the simulation to the 
experimental results of [103, 104|| (see Figs. 3 of [103,104]). Here, i?gi|(film) denotes the average 
of the parallel component of R g (z C m) over the whole film and R g n = 2i? g /3 is the corresponding 
bulk value. The ratio R s \\ (film)/i? g ii , plotted versus D/R g , indicates that the average behavior of 
the whole film becomes bulk- like if D > 5R g , as also observed in |103[| . The dotted horizontal line 
shows the statistical uncertainty of Rg\\ (horizontal dashed line). However, the overall expansion 
with decreasing D/R g is weaker in the simulation than in |103(| . 



FIG. 5. The main figure shows the variation of the chain-density profile, p c (z cm ), with the 
distance, z cm , of the chain's center of mass from the (left) wall in the high-temperature liquid 
state of the films (T = oo, i.e., only excluded volume interactions are effective). Since the profile 
is symmetric around the middle of the film, only one half is shown. Different film thicknesses, 
ranging from D = 6 (« 1.5R g ) to D = 60 (f=a 15R g ), are compared (R g (T = oo) ~ 3.68 = bulk 
radius of gyration; vertical dashed line). The small horizontal circles (•) indicate the bulk value 
Pc.bulk = 6.6 x 10~ 3 (= K/L 2 D; see Sect. ||). The inset illustrates the temperature dependence 
of Pc( •^cm) for D — 6. Here, the profile over the whole film is shown. Lines are guides to the eye 
only. The temperature varies from the high-T, liquid state (T = oo) to the supercooled state of 
the melt (T = 0.18). All lengths are measured in units of the lattice constant, and temperature is 
measured in units of e/ks (see Sect. [H]and Fig. [I]). 

FIG. 6. Temperature dependence of the monomer density profile p m (z) for D = 30 ~ 7.5R g . 
z denotes the distance of the (lower left corner of the) monomer from the (left) wall. The bulk 
radius of gyration R g increases from R g ~ 3.68 at T = oo to R g ~ 4.64 at T = 0.2 in the studied 
temperature interval. T = oo and T = 0.2 are characteristic temperatures of the high-temperature, 
liquid state and the supercooled state of the melt, respectively. The horizontal circles (•) indicate 
the bulk value /O mi buik = 6.6 x 10~ 2 (= NK/L 2 D; see Sect. |H|). Lines are guides to the eye only. 
All lengths are measured in units of the lattice constant, and temperature is measured in units of 
e/k B (see Sect. and Fig. |l|). 

FIG. 7. Time dependence of the incoherent scattering function <&q(t) for various temperatures 
ranging from the normal liquid state (T = 0.35) to the supercooled state (T = 0.2) of the melt. 
$q(t) is shown at the maximum of the static structure factor (i.e., at q = 2.94 corresponding to 
a distance of 2 (lattice constants) in real space ~ length scale of a bond). Two film thicknesses, 
D = 6 (w 1.5i? g ; depicted by symbols) and D = 30 7.5R g ; depicted by lines) are compared. 
The correspondence between temperature and line type is indicated in the figure. For the symbols 
it is: o (T = 0.35), * (T = 0.26), A (T = 0.23), □ (T = 0.2). Time is measured in Monte Carlo 
steps [MCS], i.e., the time needed to give all monomers the chance to move once. 
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FIG. 8. Incoherent intermediate scattering function $g(i) versus scaled time t/r q . The scaling 
time is defined by the condition: ^ s {r q ) = 0.629. The g-value corresponds to the maximum of the 
static structure factor (q = 2.94), as in Fig. |^. Three different temperatures are shown (T = 0.35 = 
normal liquid state of the melt; T = 0.23,0.2 = supercooled state). The symbols represent the 
data of the thin film D = 6 (~ 1.5i? g ) and the lines those of the thick film D = 30 (~ 7.5R g ). The 
correspondence between temperature and line type is indicated in the figure. For the symbols it is: 
o (T = 0.35), * (T = 0.23), □ (T = 0.2). A simple exponential, exp(-0.46i/r|), is also shown to 
illustrate that the decay of $^(i) is always stretched. The factor a = 0.46 in the exponent results 
from the requirement $g(r|) = exp(— o) = 0.629. 

FIG. 9. Temperature dependence of the relaxation times of $g(i) for two film thicknesses, 
D = 6 (~ 1.5-Rg) and D = 30 (~ 7.5R g ), and the bulk (= unconfined system of linear dimension 
L = 60). The relaxation time is defined by $q(Tg) = 0.629 (see Fig. |8|). The lines represent fits 
to the Vogel-Fulcher equation [see Eq. (||)]. The prefactor A q of Eq. (|5|) should be proportional 
to Tq(T = oo). Since the variation of r q (T = oo) with D is weak compared to that with T, it is 
possible to fit all curves with a constant prefactor At = 1.03 db 0.03. Furthermore, the activation 
energy can also be taken as constant: E q = 0.855 ± 0.005. Only To depends on D: T$ = 0.119 
(D = 6), T = 0.128 {D = 30), T = 0.136 (bulk) (error bars for all D: T + 0.001/- 0.002). The 
bulk results are compatible with earlier findings |123[ |. 

FIG. 10. Time dependence of the auto-correlation function of the fifth Rouse mode, $ 5 5 ||(i), 
measured parallel to the wall. Different temperatures are shown, which range from the normal 
liquid state (T = 0.35) to the supercooled state (T = 0.2) of the melt. The correspondence 
between temperature and line type is indicated in the figure. For the symbols it is: o (T = 0.35), 
* (T = 0.26), A (T = 0.23), □ (T = 0.2). Two film thicknesses, D = 6 (« 1.5i? g ; depicted by 
symbols) and D = 30 (~ 7.5R g ; depicted by lines) are compared with each other and with the bulk 
(indicated by x for T = 0.35 and T = 0.2). Time is measured in Monte Carlo steps [MCS]. 

FIG. 11. Same data for $55 »(t) as in Fig. |l| but plotted versus scaled time t/r^. The scaling 
time is defined by the condition: ^ss.^ts) = 0.629. The symbols represent the data of the thin 
film D = 6 (~ 1.5i?g) and the lines those of the thick film D = 30 (~ 7.5i? g ). The correspondence 
between temperature and line type is indicated in the figure. For the symbols it is: o (T = 0.35), 
A (T = 0.23), □ (T = 0.2). Furthermore, the results for the bulk are also included for T = 0.2 (x) 
and T = 0.35 (•). A simple exponential, exp(— 0.46£/t5), is also shown to illustrate that the decay 
of $55,11 (t) is always stretched. The factor a = 0.46 in the exponent results from the requirement 
$55,11(75) =exp(-a) = 0.629. 

FIG. 12. Temperature dependence of the relaxation times, T5(T), of the fifth Rouse mode 
$55,11 (t). The relaxation time is defined by $55,11 (75) = 0.629. The ordinate is divided by the infinite 
temperature result T5(oo) to illustrate the influence of film thickness with decreasing temperature. 
Two film thicknesses, D = 6 (~ 1.5R g ) and D = 30 (~ 7.5R g ), are compared with the bulk (= 
unconfined system of linear dimension L = 60). 
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FIG. 13. Schematic sketch of the influence of film thickness and temperature on static proper- 
ties. A good approximation of the instantaneous structure of a chain is an ellipsoid. The chains 
are therefore represented by ellipsoids in the figure. Adjacent to a wall the two largest principal 
axes of the ellipsoid are aligned parallel to the wall, whereas the smallest axis points away from 
it in perpendicular direction. Furtherfore, the chain is distorted (i.e., flattened). Both effects are 
entropically unfavorable so that the number of chains close to the wall is small at high temperature. 
If the film is sufficiently large, there is a bulk-like inner region, where the chains can freely reorient 
(indicated by dashed ellipsoids in the left sketch). With decreasing film thickness the bulk-like 
region shrinks more and more, leading finally to parallel oriented chains also in the middle of the 
film if the thickness is D < 4i? g . At high temperature the walls are still avoided in such a thin film 
so that a high concentration of chains occurs in the middle (indicated by the three ellipsoids in the 
right sketch instead of two in the left). Furthermore, the chains stiffen at low temperatures. This 
additionally favors orientation parallel to the walls, leading to a sequence of highly and weakly 
populated lattice layer for our model. 

FIG. 14. Illustration of a possible influence of confinement on the dynamics of a supercooled 
liquid. Two types of confinement are compared with the bulk: completely smooth walls and rough 
walls which are perfectly adapted to the liquid structure. In both cases, we assume that the walls 
are neutral. There is no strong attraction which leads to adsorption of the particles. Let us consider 
the bulk first and imagine that the particle O moves from its dashed position in direction of the 
large arrow. It opens space which its neighbors try to occupy. The dark-shaded particle has to 
compete with all other neighbors which are simultaneously moving in the space formerly occupied 
by O- This exerts a friction on the dark particle, which is partly absent if it is situated next to a 
completely smooth wall (panel in the middle). On the other hand, the friction should be strongly 
enhanced if the particle is next to a rough wall exhibiting cavities, into which it perfectly fits (right 
panel). This contraint does not relax away contrary to the bulk. So, one can expect the mobility of 
dark particle to decrease with respect to the bulk for such a rough wall (see p4] , f4lf | , for instance), 
but to increase for a completely smooth wall which acts like a lubricant compared to the bulk. 
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